Exact distribution of the maximal height of p vicious walkers 
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Using path integral techniques, we compute exactly the distribution of the maximal height Hp of 
p nonintersecting Brownian walkers over a unit time interval in one dimension, both for excursions 
(p- watermelons with a wall) and bridges (p- watermelons without a wall), for all integer p > 1. For 
large p, we show that (Hp) ~ y/2p (excursions) whereas (-ffp) ~ ^/p (bridges). Our exact results 
prove that previous numerical experiments only measured the pre-asymptotic behaviors and not the 
correct asymptotic ones. In addition, our method establishes a physical connection between vicious 
walkers and random matrix theory. 
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Introduction. Since the pioneering work of de 
Gennes [l|, followed up by Fisher [2], the subject of vi- 
cious (non-intersecting) random walkers has attracted a 
lot of interest among physicists. It has been studied in 
the context of wetting and melting (2| , networks of poly- 
mers [3| and fibrous structures ll|] , persistence proper- 
ties in nonequilibrium systems [4| and stochastic growth 
models [5|, |6j . There also exist connections between the 
vicious walkerproblem and the random matrix theory 
(RMT) lZI, Ij, p , including for instance Dyson's Brown- 
ian motion [10|. These connections to RMT have rekin- 
dled recent interest in the vicious walker problem and 
have led to new interesting questions. However, despite 
extensive recent mathematical literature on the subject, 
the connections to RMT have so far been established us- 
ing mostly combinatorial approaches. Given the non- 
intersection constraint in the vicious walker problem, it 
is natural to expect a free Fermion approach to make its 
connection to RMT physically more explicit. The aim of 
this Letter is to present such an approach which, in ad- 
dition, allows us to derive a variety of new exact results 
in the vicious walker problem. 

Physically, one dimensional vicious walkers play an 
important role in describing the elementary topological 
excitations in the p x 1 commensurate adsorbed phases 
close to the commensurate-incommensurate (C-IC) tran- 
sition Uj. In the commensurate phase the elementary 
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excitations are pairs of dislocations at a given distance 
with p nonintersecting domain walls emerging from one 
and terminating at the other. This is just a 'watermelon' 
configuration of p nonintersecting Brownian bridges (see 
Fig. lb). The sizes of such defects and their fluctu- 
ations become important near the phase transition. An 
important quantity that characterises the transverse fluc- 
tuations of the defect is the maximal height of the p vi- 
cious walkers in a fixed time (here time signifies the fixed 
longitudinal distance between the pair). Such extreme 
value questions have recently been studied extensively 




FIG. 1: a) : 4-watermelons with a wall, b) : 4-watermelons 
without wall, c) : illustration of the method to compute 
F4,{M) using path integral techniques where appropriate cut- 
offs ei's have been introduced. 



for a .single Brownian bridge or an excursion (with cer- 
tain constraints) in the context of the maximal height 
of a fluctuating interface 1^, |l3[. In this Letter, we ob- 
tain exactly the distribution of the maximal height for p 
nonintersecting Brownian bridges and excursions. 



Motivated by the geometry of elementary excitations 
discussed above, we thus focus on "watermelons" config- 
urations (see Fig. [T]a) and b)) where p non-intersecting 
Brownian walkers xi(t) < • • • < Xp(r), starting at at 
time T — 0, arrive at the same position at r = 1. We 
consider both "p- watermelons with a wall" (Fig. [T]a)), 
where the walkers stay positive in the time interval [0, 1] 
and "p- watermelons without wall" (Fig. [l]b)) where the 
walkers are free to cross the origin in between. Our main 
focus is on Hp, the maximal height of the top walker in 
[0, 1], Hp = Max^[a;p(T), < r < 1]. 

In particular, we are interested in the cumulative dis- 
tribution Fp (M) = Proba. [Hp < M] and in the moments 
(Hp). For p = 1, there exist well known results[lj|, e.g. 
(Hi) — ^7r/2 for an excursion, or (Hi) — \/tt/S for a 
bridge. Recently, Bonichon and Mosbah (BM) [l5|, us- 
ing an algorithm based on exact enumerative formulas 



|16l |. conjectured, from numerical simulations, that for 
p > 1, (Hp) aum — -\/1.67p — 0.06 for watermelons with a 
wall and (ffp)num — \/0-82p— 0.46 for watermelons with- 
out wall. These results stimulated several recent works 



171 . llSl . Iig . l2Ci [21 1 aiming at an analytical derivation of 
these estimates. 

On the other hand, exploiting the recent connection 
between watermelons and the Airy processes [^, |2^ , set- 
ting Xp{t) — Xp{t)/ y^T{l — r), one expects that, in the 
limit p ^ oo, Xp{t) = A^+p^^^^£, where A — 2^/^ (ex- 
cursions) and A — 2 (bridges), where ^ is the Airy2 pro- 
cess [9|,|2^ of a suitably rescaled time parameter. Thus in 
the large p limit, the top curve approaches a limit shape, 
Xp{t) — !■ A^/py^T{l — t). Since the maximum of the top 
curve occurs at the midpoint t = 1/2, one expects that 
for p > 1, ( ffp ) - ( Xp{T = i)) - v^ for ex- 
cursions and, similarly, {Hp) ~ ^ for bridges. These 
exact asymptotic estimates differ considerably from the 
numerical estimates of BM suggesting that the latter only 
describe the preasymptotic behavior of {Hp). However, it 
calls for an explanation why this preasymptotic behavior 
as measured by BM should be about y/1.67p and -^0.82^. 

In this Letter, we present a method based on path inte- 
grals associated to corresponding free Fermions models to 
compute exactly Fp{M). Our exact formula is useful for a 
number of reasons. It provides the exact asymptotic tails 
of the distribution of Hp which were not known before. 
For the average height, our formula explains the afore- 
mentioned discrepancy between the estimates of BM and 
the exact asymptotic behaviors of (Hp) . We show that for 
moderate values of p (preasymptotic behavior), one ob- 



-^1. 64493 • • -p for excursions and 
-^0. 822467- • -p for bridges, in nice 
agreement with BM's estimates. Finally, we show how 
our method allows for a physical derivation of the con- 
nection between p- watermelons configurations and RMT. 
Method. To calculate the cumulative distribution 
Fp{M), we use a path integral method which needs to 
be suitably adapted to this problem. Indeed one notices 
that the p-watermelons configurations described above 
(see e.g. Fig. [T]a) and b)) are ill defined for systems in 
continuous space and time. For such Brownian walks, it 
is well known that if two walkers cross each other once, 
they will re-cross each other infinitely many times imme- 
diately after the first crossing. Therefore, it is impossible 
to enforce the constraint Xi{0) = a;i+i(0) = and si- 
multaneously forcing Xi{T) < a;j+i(T) immediately after. 
The cleanest way to circumvent this problem is to con- 
sider discrete time random walks moving on a discrete 
one-dimensional lattice (so called Dyck path) : this is 

l3, 111, Si- By taking 



tains {Hp) ex n 
{Hp) ex TTv^VlS 



the method used in Ref. [l6|, 
the diffusion continuum limit, one would then arrive at 
non intersecting Brownian motions [23|. This method is 
however mathematically cumbersome. Alternatively, fol- 
lowing Ref. [13, |2J| , we can go around this problem by 
assuming that the starting and finishing positions of the 



p walkers are < ei < ... < e^ (see Fig. [T]c)). Only at 
the end we take the limit ei ^ and show that it is well 
defined. In addition, in order to compute Fp{M), we put 
an absorbing hard wall at M such that 



Fp{M) = lim 



N{e, M) 



N{e, M -> 00) 



(1) 



where e = ei, - - - , Sp and N{e, M) is the probability that 
the p Brownian paths starting at < ei < ... < ep at 
T — Q come back to the same points at r = 1 without 
crossing each other and staying within the interval [/i, M], 
with p, = Q for excursions and p, — > — cx) for bridges. This 
procedure is depicted in Fig. [T]c). 

The probability measure associated to p unconstrained 
Brownian paths xi{t)^ ..,Xp{T) over the time interval 

[0, 1] is proportional to exp [—i ^^^^ Jg (^) '^''']- Here, 
we have to incorporate the constraint that they stay in 
the interval \p,M]. Therefore one can use path-integral 
techniques to write iV(e, M) in Eq. H]) as the propagator 



N{e,M) 



-H^ 



(2) 



with Hj^i = 'J2^i=i[^'i^~^^i^i)]j where V{x) is a confin- 
ing potential with V{x) — ii x £ [/i, M] and V{x) — 00 
outside this interval. Denoting by E the eigenvalues of 
Hm and \E) the corresponding eigenvectors one has 



N{e,M)=Y,\'^Eie)\ 



^e-^ 



(3) 



where we introduced the notation (x|£') — ^£;(x). Im- 
portantly, to take into account the fact that we are con- 
sidering here non-intersecting Brownian paths, the many 
body wave function ^b(x) = "^^ixi, ..,Xp) must be 
Fermionic, i.e. it vanishes if any of the two coordinates 
are equal. This many-body antisymmetric wave function 
is thus constructed from the one-body eigenfunctions of 
Hm by forming the associated Slater determinant. 

Watermelons with a wall. In that case /i = 
and the one-body eigenfunctions are given by (j)n{x) ~ 



2 „;„ mrx 
M ''^^ M 

Therefore one has 



sm -^^^ with discrete eigenvalues ^^^ , n e 



2M2 ' 



-= det 0„.(ej) , E 



*£;(e) = ^ det , 



2Af2 



n 



(4) 



where we use the notation n^ — J^^^i^h "» '= '^*- 
From this expression (|4]), one checks that, in the limit 
£1 —f 0, - - - ,ep — *■ 0, powers of e^'s cancel between the 
numerator and the denominator in Eq. ([T]), yielding 



FpiM) 



A^ri 



]\,f2p^^p 



E [-(") 



(5) 



5(n)= n K'-"')!!"^' 
i<j<fc<p j=i 



where Ap, a constant independent of M, is determined instance, F2{M) = 1 — 4M^e" 



-2M^ 



-4M^ 



by requiring that lim.M^oo F{M) = 1. It can be eval- 
uated using a Selberg's integral [2^ yielding Ap = 
7rV+P/[2P'-p/2j]P-ir(2 + j)r(|+j)]. Forp = 1, our 
expression gives back the well known result for a Brow- 
nian excursion |26| . For p = 2, we have checked, using 
the Poisson summation formula that our expressions in 
Eq. ([5]) yield back the result of Ref. [18[. For generic p, 
the probability distribution function (pdf) F'{M) is bell- 
shaped, exhibiting a single mode. At variance with previ- 
ous studies 13, llSl , our expression ^ is easily amenable 
to an asymptotic analysis for small M. Indeed, when 
M -^ 0, the leading contribution to the sum in ([5]) comes 
from n.i = i and its p\ permutations, yielding for A/ — > 



Fp{M) 



_/\,f2p2+p 



^p(p+l)(2p+l) 



(6) 



where ap can be explicitly computed, yielding for in- 
stance a2 = 127r^. For large M, one can use the Poisson 
summation formula to obtain 1 — Fp{M) ex exp(— 2M^). 
From the distribution in Eq. ([5]), one can compute the 
moments of the distribution (ffi). For p > 2, one obtains 
that {Hp) can be expressed in terms of integrals involving 

the Jacobi theta function ^{u) — Yl^=-oo e"'^" " and its 

b y a simpler physical deriva- 



ivatives, tnus recovering, by a simpler pnysical deriva- 
1, the results of Ref. [llJSi, IH, [l^. In particular 



derivatives, thus recovering 
tion, 

one has {H2) = 1.82262... [17|. For moderate values of 
p, one observes that the main contribution to the aver- 
age {Hp) = /g MF'{M)dM comes from relatively small 
M where F'{M) is dominated, as before in Eq. ([S]), by 
the terms where rii = i and its p\ permutations. It is 
easy to see that the pdf, restricted to this first term ^ 
exhibits a maximum for M* ^ TTy/p/6. Therefore, one 
expects that (Hp) ^ M* = a/1.64493 • • -p, in good agree- 
ment with the estimates of BM [151]. For larger values of 
p the average {Hp) picks up contributions from larger 
values of M where F'{AI) can not be approximated by 
a single term as in Eq. ^ and therefore the estimate 
of BM ceases to be correct. Instead, one has the exact 
asymptotic behavior {Hp) ~ \/2p for p ^ 1, which can 
be obtained directly from our formula in Eq. ([5]) !29i] . 

Watermelons without wall. In the case of Brownian 
bridges, one can apply the same formalism as above 
© - © with ^ -^ - ^, t.e. Hm = Er=i ^^■ 
In that case, the one-body eigenfunctions are given by 
tpf^[x) — \ — sin [k{M ~ x)] with a continuous spectrum 

Ek = fc^/2, k e R+. Therefore, *_B{e) entering the ex- 
pression of 7V(e, M) in Eq. ([3]) is formally given by Eq. @ 
where (^„. is replaced by ipk. and E = ^. One obtains 



Fp{M) 

0p(y) = 



Bp 



dyi 



dVp 



"2M2 



0p(y)' 



det y\ 



•■' ^cos(2/,+jf), (7) 

where Bp = 2"^^ /[{2ti)'p''^ Jlj^i r(j + 1)]. This yields, for 



From Eq. 
as 



(I?!), one obtains the asymptotic behavior for A4 

Fp{M) ex KP^+P , (8) 

whereas for large M one has 1 — Fp{M) (x exp (— 2M^). 
As in the case of watermelons with a wall, the pdf F'{M) 
is also bell-shaped with a single mode. Notice however 
that the presence of the wall has drastic effects on the 
small M behavior of Fp{M) (see Eq. dH) and Eq. §^) 
whereas, as expected, it has less influence for large M. 
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FIG. 2; a) : Contour plot of 02(2/1,5/2) = 2/2 sinyl cos j/2 — 
yi cos j/i sin j/2 given in ([7]). It exhibits saddles for (j/i,j/2) = 
(±7r/2, ±7r) and symmetric points obtained by permutations, 
b) Plot of {Hp) as a function of p. The dotted line is the 
estimate from BM ]l3^ . The quality of this estimate for p ^ 10 
has its origin in the saddles of Op(y) shown, for p = 2, on the 
left panel. For larger values of p one has instead {Hp)'^ ex p. 

From Fp{M) in ([7]), one computes the moments {Hp), 

yielding {H,) = l±^^ or (i/3) = ^^±^^#-^7^, 
recovering (to leading order) recent results obtained by 
rather involved combinatorial techniques [2C| . 

To make contact with BM's estimates, one first focuses 
on p = 2 and notices that 62(2/1, 2/2) in Eq. ([7|) exhibits 
saddles for yi = ±7r/2, j/2 = iT^ and for symmetric points 
obtained by permutations: this is shown in Fig. [5] a) . In 
fact this property can be generalized to higher values of 
p and one can show that 0p(y) has saddles which are lo- 
cated around yi — ±7r/2, j/2 = iTr, ■ • ■ , 2/p = ±p7r/2 and 
the points obtained by permutations. Of course Op(y) 
develops saddles for higher values of y^ but their weights 
are exponentially suppressed in Eq. ([T]). For moder- 
ate values of p, one expects that {Hp) is dominated by 
these saddles j/i = ±7r/2,?/2 — ±7r, • • • ,?/p = ±p7r/2. 
Therefore performing a saddle point calculation, one has 

Fp{M) (X e P ^\^) , with x{y) = log?/ + 7r2/(24y2), 
which has a minimum for y* = 7r/vT2. This yields 
{Hp) ^ TT^Jp/12 = 1/0-822467 •••p, in good agreement 
with the estimates of BM [15J . For larger values of p one 
expects that {Hp) picks up contributions from larger val- 
ues of M where Fp{M) can not be reduced to these first 
saddles. In Fig. [2]b), one shows a comparison between 
the exact value of {Hp)'^ computed from Eq. ([7]) and the 
estimate of BM. This clearly shows that the estimate of 
BM correspond to the pre-asymptotic behavior. Instead, 
for large p, one expects here {Hp) ex y/p. 



Extension of the method. The method presented here 
can be used to derive many other results. As an interest- 
ing example, showing explicitly the connection between 
watermelons and RMT, we compute the joint probability 
distribution Pjoint(a^i,' • ■ ,Xp,T), first for p bridges. Fol- 
lowing the same steps as above, Eq. ((H)-©, and using 
the Markov property of Brownian paths, one has 



asymptotic behaviors for the average (Hp). Besides, we 
obtained a quantitative description of the pre- asymptotic 
regime actually measured in the numerical experiments 
of BM. We hope that the path integral method presented 
here, which is rather general, and the precise connection 
to RMT will allow further future studies. 
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lim 






(9) 
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with Hq — X]r=i "J"^;^- One can show that powers of 
Ei's cancel between the numerator and the denominator 
in (O, yielding Pjoint(x, r) ex (3(x, t)Q{x, 1 - r) with 



Q{j^,t)^ f dkHih - k,)e 



■2 det e 

l<77i,n<p 



{iXTnk-n) 



,(10) 



where J dk = /_ dki • • • /_ dkp. After some algebra 
to evaluate the integrals in Eq. p^ one finally obtains, 
for p- watermelons without wall 

Pjoi„t(x,T) = Z-V(r)-P' n(^* - 2^.)'e"^ , (11) 



with aij) — \/t{\ — t) and Zp a normalization constant. 
This expression in Eq. pT|) shows that this joint proba- 
bility is exactly the one of the eigenvalues of the Gaussian 
Unitary Ensemble of random matrices (GUE) |7l.l8lllCl|. 
In particular, for p 3> 1, defining the rescaled variable 
■q — y^pi/6( "^ \" , — \/2p)j one obtains that the cumula- 
five distribution of r/ is given by Proba[77 < x\ — J-2 {x) , 
the Tracy- Widom distribution for /3 = 2 [irf . 
For excursions, a similar calculation shows that 



Pjoi„t(x,r) = Z'-V(r)-^(^''+i)[S(x)]^ 



e 2<t^(t) 



(12) 



where S(x) is defined in ([5]) and Z' a normalization con- 
stant. Hence the joint distribution of yi = xf/2cr^(T) 
is formally identical to the distribution of the eigenval- 
ues of Wishart matrices [25| with M — N — i, and 
N = p. In that case, from the results for the largest 
eigenvalue of Wishart matrices we conclude that for p ^ 
1, the cumulative distribution of the rescaled variable 
C = 22/3pi/6(_^lL _ 2^) is again given by J^2ix) Q. 
Conclusion. To conclude, using methods of many-body 
physics, where appropriate cut-offs e^'s have been intro- 
duced (see Fig. [T]c)), we have obtained exact results for 
the distribution of the maximal height for p- watermelons 
with a wall ([5|) and without wall ([7]), which is physically 
relevant to describe the geometrical properties of disloca- 
tions arising in p x 1 commensurate adsorbed phases close 
to the C-IC transition. Our expressions explain the dis- 
crepancy between the estimates of BM [15| and the true 
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